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Abstract 

We show that the study of the statistical properties of the scattering 
matrix S for quantum chaotic scattering in the presence of direct pro- 
cesses (characterized by S ^ 0, S being the average S matrix) can be 
reduced to the simpler case where direct processes are absent (S = 0) . 
Our result is verified with a numerical simulation of the two-energy 
autocorrelation for two-dimensional S matrices. It is also used to ex- 
tend Wigner's time delay distribution for one-dimensional S matrices, 
recently found for S — 0, to the case S ^ 0; this extension is verified 
numerically. As a consequence of our result, future calculations can be 
restricted to the simpler case of no direct processes. 

The problem of chaotic wave scattering is of great interest in nuclear, 
molecular and mesoscopic physics, in optics and in the microwave domain. 
Common features have been found for system sizes spanning such a large 
range because of the generality of the phenomena involved. 

In quantum-mechanical scattering problems with a chaotic classical dy- 
namics one aims at studying the statistical properties of the scattering ma- 
trix S. Recently, this study has been further motivated by experiments on 
quantum electronic transport in mesoscopic systems Q. 

The one-energy statistical distribution of the S matrix has been de- 
scribed successfully by an information-theoretic model that incorporates 
precisely the physical information that is relevant for a wide class of sys- 
tems [Q, pi |4|. That information specifies: 1) General properties: i) flux 
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conservation (giving rise to unitarity of the 5 matrix), ii) causality and the 
related analytical properties of S(E) in the complex-energy plane, and Hi) 
the presence or absence of time-reversal and spin-rotation symmetry, that 
determines the universality class: orthogonal, unitary or symplectic (desig- 
nated as f3 = 1,2,4) H ||] and restricts further the structure of 5: unitary 
symmetric, unitary or unitary self-dual, respectively. 2) A specific prop- 
erty: the ensemble average (5), identified with the energy average 5, also 
known as the optical S matrix |7j , which controls the presence of prompt, or 
direct processes in the scattering problem. In this procedure one constructs 
the statistical distribution of 5 using only the above physical information - 
expressible entirely in terms of 5 itself- without ever invoking any statistical 
assumption for the underlying Hamiltonian, that never enters the analysis. 
The resulting 5-matrix distribution, known as Poisson's kernel, reproduces 
well the statistical scattering properties of ballistic cavities with a chaotic 
classical dynamics Qj. 

The joint statistical distribution of the S matrix at two or more energies 
has escaped, so far, an analysis within the philosophy described above (some 
aspects of the two-point problem have been studied assuming an underlying 
Hamiltonian described by a Gaussian ensemble, as in Refs. [||, J9[). An 
approach coming close to that philosophy was initiated in Ref. [[1C| ], with 
a study of the simplest quantity of a two-point character: the statistical 
distribution of the time delay -that involves the energy derivative of the 5 
matrix- arising in the scattering process [11|. The physical motivation of 
Ref. [10 1 was the study of the electrochemical capacitance of a mesoscopic 
system [l^ ]. The analysis was done for a cavity attached to one lead that 
can support only one open channel (so that 5 = e tS is a 1 x 1 matrix), for 
arbitrary (3 and 5 = 0. It was based on a conjecture by Wigner |l3|], that 
assumes invariance of the statistics of the poles and residues of the related 
K matrix under the same transformation 5 — > e 1 ^/ 2 Se 1 ^/ 2 that defines the 
invariant one-point measure [|14|] which, in the one-channel case, is d9/2ir. 
That study was generalized to an arbitrary number of channels N in Refs. 



[ 15| , [L6fl , again for 5 = 0: Wigner 's conjecture was formulated and proved 
as the invariance of the fc-point probability distribution of the A^-channel 5 
matrix under the transformation that defines the invariant measure for the 
universality class /3. 

Of physical importance is the case 5 ^ 0: it corresponds to situations 
where direct processes are not negligible. For example, the case 5 = 
describes chaotic cavities with ballistic point contacts, while the coupling to 
leads containing tunnel barriers produces direct reflection and thus 5^0 
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[ I , [l?]] . A more complex combination of direct processes is described in Ref . 
If- 

In the present article we find a transformation that relates the Appoint 
distribution of the iV-channel S matrix for the case S ^ to that for 5 = 0, 
thus relating the problems in the presence and in the absence of direct pro- 
cesses. In Ref. |fE8[ ) the problem with a nondiagonal S was reduced to that 
with a diagonal and real S. Here we show that the problem can be reduced 
further to one characterized by S = 0. That transformation, used in the 



past in relation with one-point functions [19], contains only the S matrix 
and its average S. It is used here to extend to the case 5^0 the time-delay 
distribution of Ref. |l0| and is verified in this case by comparing our results 
with some known ones and with numerical simulations. It is also verified 
numerically for the two-point correlation function of two-dimensional S ma- 
trices. We believe that our result is appealing because of its conceptual 
simplicity and may open the way for further studies in this field. 

We first summarize the information-theoretic model of Refs. ||, 0]. 
The starting point is the measure d/j,^'(S) for n x n S matrices, that re- 
mains invariant under the symmetry operation for the universality class (5 
[]|, ^, 14 1 . The ensemble average of S vanishes for the invariant measure, 



and so does the prompt component. An ensemble for which (S) = S / 
contains more information than the invariant-measure: it can be con- 
structed as (S)dfi^ (S) and the information associated with it is defined 



as / p-^\S)ln p^ 1 (S) dfi^(S). We assume to be far from thresholds and 
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recall that S(E) is analytic in the upper half of the complex-energy plane 
(causality) [^0| . The study of the statistical properties of S is simplified by 
idealizing S(E), for real E, as a stationary random function satisfying the 
condition of ergodicity, which in turn implies the equivalence of energy and 



ensemble averages [21|. We thus have the analyticity- ergodicity requirements 
(AE), that imply {S^ b S q d ■ ■ •) = (S a t,) p {S C( i) q • • •, so that averages of prod- 
ucts of /S-matrix elements not involving complex conjugation can be written 
in terms of the matrix (S) = S. The one-energy probability density 

P S (D) V P I det(I-SSt) |/3n+2-/9 ' W 

where Vg is a normalization constant, is known as Poisson's kernel: it satis- 
fies the AE requirements Q and the associated information is less than or 
equal to that of any other probability density satisfying AE for the same S 
||. In Eq. (||), the elements of S are assumed to be complex numbers for 
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/3 = 1, 2 and quaternions for j3 = 4; for the definition of the determinant of 
a quaternion matrix, see Ref. [||], p. 126 . 
Now consider the transformation 

S = t' 1 - 1 (S-r 1 )(l-r{S)- 1 ti (2) 

where r\, t\, r[, t'i are the n x n blocks of the matrix 



(3) 



which has the symmetry associated with the universality class (3. One can 
prove the following statement j2|, [l(| [l?], |l9| : if the one-energy distribution of 
S is Poisson's measure (||) and we identify r\ in the transformation (|2|) with 
S, the one-energy distribution of So is the invariant measure dfi^\So). In 
other words, the transformation (|2|), with r\ = S, transforms the problem 
with direct processes to one without (So = 0), as far as the one-energy 
distribution goes. 

Now suppose that at every energy E we subject S(E) to the transforma- 
tion (|2|), always with the same S. We prove below the following statement: 
the joint statistical properties of the transformed Sq(Ei), Sq(E2) , etc. are 
precisely the ones associated with the problem without direct processes, char- 
acterized by So = 0. In other words, the above transformation relates S(E), 
understood as a stationary random function of energy, for a problem with 
direct processes to one without. 

We prove the above statement for S unitary and symmetric ((5 = 1), the 
proof for the other symmetries being similar. Consider first the case of S 
diagonal and real. S can be written in terms of the real amplitudes 7Aa's 
and the energy levels E\s, as 

S(E) = [l + tK(E)][l-iK(E)]- 1 , (4) 
where the K matrix is given by 

K ab (E) = J2f^%- (5) 
- E\ — E 

The 7Aa' s are uncorrelated Gaussian variables and the E\s follow the statis- 
tics of the Gaussian Orthogonal Ensemble || with average spacing A. 
S(E) is given by @ 



S(E) ab = (l-y aa )(l + y bb )' 1 5 ab , (6) 
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where y aa = 7r(7^ a )/A. Just as S and K are related by Eq. (||), the 
transformed Sq, Eq. (|2|), can be written in terms of a matrix Kq given by 

K Q = i(l-So)(l + S )- 1 . (7) 

Substituting the transformation @ into (0), with S given by we finally 
find 

(K )a6 = Vaa /2 K ab y; b 1/2 . (8) 

— 1/2 

Thus, the 7Aa's are rescaled by a constant (7a q — > yaa 7Aa), which is the 
appropriate one to ensure Sq(E) = 0. Therefore, under the transformation 
@, the 7Aa's retain their Gaussian distribution, while the energy levels E\s 
are unchanged. We conclude that the statistical properties of the trans- 
formed random function Sq(E) are the same as those of a statistical S(E) 
matrix without direct processes. This proves our statement for S diagonal 
and real. The case of arbitrary S can now be reduced to the above one using 
the results of Ref. @. 

As a first application of the above statement, we extend to S ^ the 
time-delay distribution found in Ref. |T(| for one spatial channel, S = 
and arbitrary (3. We consider Eq. (|2|) for n = 1, with r\ = S; omitting the 
unimportant phase factor t*/t\, we have 

S (S) = (5-5)(l - S*S)-\ (9) 

where So and S can be written as Sq = exp(i6>o), 5 = exp(i#). 

We are interested in the time delay JTB] 0' = d9/dE; we express it in 
terms of 6q and 0' = dO^/dE as 



0' = dO/dE 



1 



_ 2 

5 



1 + S*e* 



dSo/dS = /(^o)^o- (10) 



To find the distribution of 9' we need the joint distribution of 0q and 9' , 
Pq(9q, 9q), which, being the one for no direct processes, factorizes as ]l5| 

vf\9 ,9' )=vi ) (e' )/2n. (11) 



It is convenient to express Eqs. (|10| ) and ( |ll[) in terms of the variables u, uo 
defined as 

u = 2tt/9'A, u = 2tt/9' A , (12) 
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where A and Ao denote the average level spacing for the problems described 
by S and So, respectively. The above proof shows that A = Ao, since the 
energy levels E\ are unchanged by the transformation (^). Thus, using ( p"2| ) 



and no = 2tt/9q'A, Eq. (10) becomes 

u = u /f(6 ). (13) 
The joint distribution of the statistically independent variables 6q,uq [see 



Eq. (11)], needed to find the distribution of u, is given by 

P o (/3) (0o,n o ) = P o (/3) K)/2vr, (14) 

where Pq^\uq) was obtained in Eqs. (16), (17) of Ref. [[[(]] [and denoted 
there by P(u)] as 

F ? W = fl[ a » /V(M " ' (15) 

From @, (0), (|l|) we write the distribution of u, (u) = (6[u- u //(6>o)]), 
as 

if \u) = ±.f f(e )P^(f(o )u) d e = ^^/ ,2 fJ lf(e )] l+m e-^d9 . 

For the dimensionless time delay t = 1/u = 6'A/2ir we finally find the 
probability density 



(f) 



P/2 

2tt 



~ % / * e-frWdBo. (17) 

2vrr (|) r 2+ 2 Jo 



The calculation of w-jt (r) for arbitrary S and /? is thus reduced to 



} i_P) 
g - ■ 

quadratures [f(@o) is known, Eq. (|i~0|)1 and the result coincides with our 
previous one |l(| for S = 0. 

Eq. ( |T7| ) is compared in Fig. 1 (a) with a numerical simulation that 
generates an ensemble of S"s from resonances sampled from an unfolded 
Gaussian ensemble for (3 = 1, the coupling amplitudes to the channel being 
independent Gaussian variables whose variance ensures S = 1/2. The agree- 
ment is very good. Such a comparison was also made for /3 = 2, 4, although 
it is not illustrated here. For j3 = 2, Ref. || gives w^\t) in analytical form 
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in terms of Bessel functions: we could not prove analytically the equivalence 
with our result; however, numerically the two are indistinguishable. 

As a further verification of our statement, consider the autocorrelation 
function Gg(E) = (Sf^SniE))^ - (^(O))^ (S n (E))-g. The quantity 

i 1 2 

|c>g(i£)| was calculated by generating numerically ensembles of 2 x 2 5 
matrices, with 5 = and (1/2)1, I being the unit matrix. To the data for 
5 = (1/2)/ the transformation (|2|) was applied: the results shown in Fig. 1 
(b) are seen to be consistent with \cg =0 (E)\ 2 . 

Summarizing, we have found a transformation that relates the 5 matrix 
S(E), understood as a stationary random function of energy, for a problem 
with direct processes (5 ^ 0) to one without (5 = 0). An application was 
made to extend Wigner's time-delay distribution for one channel 5-matrices, 
from 5 = to 5 ^ 0. A number of numerical simulations was made as a 
verification of our transformation. Our result implies that future work on 
the statistical properties of the 5 matrix can be restricted to the simpler case 
5 = 0, and extended to the case 5^0 —corresponding to more complex 
scattering systems with direct processes— using the procedure described in 
this paper. 
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Figure 1: (a) The distribution of time delays for one channel and in the 
presence of direct processes (5 = 1/2), for = 1. The crosses are propor- 
tional to the theoretical probability density of Eq. ( |l7|) , that was integrated 
numerically. The points with the finite-sample error bar are the results of 
the numerical simulation described in the text. The agreement is excellent, 
(b) The square of the autocorrelation function of the 5n element of a two- 
channel 5 matrix for (5 = 2 as a function of the energy separation, obtained 
from a numerical simulation. The open circles and diamonds correspond to 
5 = and (1/2)/, respectively. The squares are the result of applying the 
transformation (||) to the data for 5 = (1/2)1: they are seen to be consistent 
with the correlation for 5 = 0. 
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